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u ■ Abstract 

We classify (l+3)-dimensional Pauli equations for a spin-i particle interact- 
ing with the electro-magnetic field, that are solvable by the method of separation 
of variables. As a result, we obtain the eleven classes of vector-potentials of the 
electro-magnetic field A(t, x) = (Ao(t,x), A(t,x)) providing separability of the cor- 
responding Pauli equations. It is established, in particular, that the necessary con- 
dition for the Pauli equation to be separable into second-order matrix ordinary 
differential equations is its equivalence to the system of two uncoupled Schrodinger 
■ equations. In addition, the magnetic field has to be independent of spatial variables. 

We prove that coordinate systems and the vector-potentials of the electro-magnetic 
field providing the separability of the corresponding Pauli equations coincide with 
Q-T those for the Schrodinger equations. Furthermore, an efficient algorithm for con- 

structing all coordinate systems providing the separability of Pauli equation with a 
fixed vector-potential of the electro-magnetic field is developed. Finally, we describe 
all vector-potentials A(t, x) that (a) provide the separability of Pauli equation, (b) 
satisfy vacuum Maxwell equations without currents, and (c) describe non-zero mag- 
netic field. 
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1 Introduction 

A quantum mechanical system consisting of a spin— ~ charged particle, moving with mo- 
mentum p in a time-dependent electro-magnetic field with the four-component vector- 
potential (Aq,A), is described in a non-relativistic approximation by the Pauli equation 
(see, e.g., gj) 

(po-eA {t,x) - [p-eA^x)^ + eaHj ip{t,x) = 0. (1) 

Here ?/>(£, x) is the two-component wave function in three space dimensions x = (x±, X2, 23), 
H = rot A is the magnetic field, and a = (<7i, cr 2 , 03) is a vector consisting of three Pauli 
matrices 
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Hereafter we use the notations 

d 

Po = i gl^ p= -iV, a = 1,2, 3, (3) 

and summation over the repeated Latin indices from 1 to 3 is implied. 

As the Pauli equation has variable coefficients, we cannot apply the standard Fourier 
transformation. The only regular way for solving (|l|) is the classical method of separation 
of variables in curvilinear coordinate systems. In this respect, a natural question arises, 
which equations of the form (|l|) are separable, namely, which potentials Aq, A allow for 
separability of the Pauli equation in some curvilinear coordinate system? One of the 
principal objectives of the present article is to provide an efficient algorithm for answering 
these kinds of questions for systems of partial differential equations. It is essentially 
based on the results of the paper 0, where the problem of separation of variables in the 
Schrodinger equation has been analyzed in detail. As the Pauli equation (p]) differs from 
the Schrodinger equation by the term eaH only, it is natural to attempt modifying the 
technique employed in in order to make it applicable to system of partial differential 
equations (PDEs) (0). 

Integrable Hamiltonian systems with velocity-dependent potentials have been studied 
for the case n — 2, i.e., in a Euclidean plane by Winternitz with co-authors ffl. Recently 
Benenti with co-authors studied the problem of separation of variables in the stationary 
Hamilton- Jacobi equation with vector-potential from a geometrical point of view. 

The problem of separation of variables for linear systems of first-order partial differen- 
tial equations such as the Dirac equation has been repeatedly addressed by Shapovalov 
and Bagrov with co-authors |6], [7j] and by Kalnins and Miller with co-authors |8|, |9], [Hj . 
They developed a symmetry approach to the separation of variables in the Dirac equation 
where separability is characterized by the existence of a complete set of first-order matrix 
symmetry operators. 

Symmetry and supersymmetry properties of the Pauli equations are studied in [O 



121 [13|. Let us also mention the paper [14|, where physical aspects of the problem of 



separation of variables in some (l+3)-dimensional Pauli equations with time dependent 
potentials are studied, and the monograph |7|], where some classes of exact solutions of 
the Pauli equation are presented. 

With all the variety of approaches to separation of variables in PDEs one can notice 
the three generic principles, namely, 

a) Representation of a solution to be found in a separated (factorized) form via several 
functions of one variable. 

b) Requirement that the above mentioned functions of one variable should satisfy some 
ordinary differential equations. 

c) Dependence of so found solution on several arbitrary (continuous or discrete) pa- 
rameters, called spectral parameters or separation constants. 

By a proper formalizing of the above features we have formulated in an algorithm for 
variable separation in the Schrodinger equation with vector-potential. Below we generalize 
this algorithm for the case of system of PDEs (|). 
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To have a right to talk about description of all potentials and all coordinate systems 
enabling us to separate the Pauli equation, one needs to provide a rigorous definition of 
separation of variables. The definition we intend to use is based on ideas contained in the 
paper by Koornwinder [ T5\ . 



Let us introduce a new coordinate system t, u a = u a (t,x), a = 1,2,3, where u a are 
real- valued functions, functionally independent with respect to the spatial variables X\, 
x 2 , x 3 , i.e.: 



det 



* 0. (4) 

a,b=l 



dx b 

For a solution to be found we adopt the following separation Ansatz: 
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ip(t, x) = Q(t, x)(p (t) Yl Va (u a (t, x), \ \ x, (5) 



a=l 

where Q, tp^, (/i = 0, 1, 2, 3) are non-singular 2 x 2-matrix functions of the indicated 
variables and \ is an arbitrary two-component constant column. What is more, the 
condition of commutativity of the matrices (p^ is imposed, namely, 

[<Pn, fv] = Witfv - VvWiL = 0, n, v = 0,1,2,3. (6) 

Note that the restriction (|6]) is an extra requirement, which narrows the class of separable 
Pauli equations. However, without this condition an efficient handling of the Ansatze of 
the form (|6|) seems to be impossible. At least, in all papers devoted to variable separation 
in a systems of PDEs the condition of commutativity is imposed (explicitly or implicitly). 

Definition 1 We say that the Pauli equation ([Jj admits separation of variables in a 
coordinate system t, uj a = u a {t,x), a = 1,2,3, if there are non-singular 2 x 2-matrix 
function Q(t,x) and four matrix ordinary differential equations 

iipo = - (Poo(t) + Pob(t)\b) <Po, 

(7) 

<Pa = (-Pao(^a) + P a b(^a)K) Pa, a =1,2, 3, 

jointly depending in an analytic way on three independent complex parameters Ai, A2, A3 
(separation constants), such that, for each triplet (Ai, A2, A3) and for each set of solutions 
<fo(t), y?i(^i) ; ^2(^2); ^3(^3) °f iW' f unc ti° n (Hi under condition ^) is a solution of (fjj. 

In the above formulas P^, v = 0, 1, 2, 3 are some complex 2 x 2-matrix functions of 
the indicated variables. 

Definition 2 Three complex parameters Ai,A2,A3 in (Qj are called independent, if the 
equality 

rank ||P ua || 3 n 3 , = 6. (8) 

II ^"11^=0 a=l \ / 



holds, whenever (po(t){pi(u}i)ip2(u>2) l P3{u3) ^0. 
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Condition (|3|) secures essential dependence of a solution with separated variables on 
the separation constants A. 

Note, that putting Q = I,u a = x a ,a = 1,2,3 in (El) yields the standard separation of 
variables in the Cartesian coordinate system. Next, choosing the spherical coordinates as 
u>i, u 2 ,uJ3 we arrive at the variable separation in the spherical coordinate system and so 
on. The principal task is describing all possible forms of the functions Q,u a ,a = 1,2,3, 
that provide separability of the Pauli equation in the sense of the definition given above. 
Solution of this problem, in its turn, requires describing the functions A , . . . , A 3 that 
enable variable separation in the Pauli equation in the corresponding coordinate system. 
More precisely, we will need to solve the two mutually connected principal problems: 

— to describe all cases of coefficients, for which the corresponding Pauli equation ([[]) 
is separable (in the sense of definition 1) in at least one coordinate system; 

— to construct all coordinate systems that allow for separation of variables (in the 
sense of definition 1) in the Pauli equation ([![) with some fixed vector-potential 

Note, that formulas ©-(|) form the input data of the method. We can change these 
conditions and thereby modify the definition of separation of variables. For instance, we 
can change the order of the reduced equations ([?]) or the number of essential parameters 
A a (a more detailed analysis of this problem for the Schrodinger equation can be found 
in [16|]). So, our claim of obtaining the complete description of vector-potentials and 
coordinate systems providing separation of variables in ([I]) makes sense only within the 
framework of definition 1. If one uses a more general definition, it might be possible to 
construct new coordinate systems and vector-potentials providing separability of equation 
(H|). But all solutions of the Pauli equation with separated variables known to us fit into 
the above suggested scheme. 

Tr ansf or mat ions 

A a ^A' a = A a (Ai,A 2 ,A 3 ), a = 1,2, 3 (9) 

under condition 

^ 0. (10) 

a,b=l 

preserve the form of relations (§)-©■ So we can regard the corresponding spectral pa- 
rameters A and A' as equivalent ones. Within the framework of this equivalence relation 
we can choose A in such a way that all matrices P^, fi, v — 0, 1, 2, 3 in reduced equations 
(0) are Hermitian ones and parameters Ai, A 2 , A3 are real numbers. 

Next, we introduce an equivalence relation on the set of all vector-potentials A (t,x), 
A(t,x) providing separability of equation ([I]), on the sets of solutions with separated 
variables and corresponding coordinate systems. 

Definition 3 We say that two vector-potentials A(t,x) and A'(t,x) are equivalent if they 
are transformed one into another by the gauge transformation 

- - df 
A ^ A ' = A + V f, Ao ^A' = A -^L, (11) 

where f — f(t, x) is an arbitrary smooth function. 
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For the Pauli equations to be invariant with respect to the above transformation, the 
wave function i/j(t,x) is to be transformed according to the rule 

iP _> $ = ip exp(ief) (12) 

Indeed, if the transformations (|TT|)-(|T2"|) in the Pauli equation ([]]) are performed one after 
another, we obtain the initial equation, provided we replace the functions A,Ao,ip with 
■V. %■>■>. 

Note that the system of PDEs (|l|) admits a wider equivalence group from the point of 
view of the standard theory of partial differential equations (Shapovalov and Sukhomlin 
fT7fl). However, this group cannot be regarded as an equivalence group within the context 
of quantum mechanics, since allowed transformations of the wave function must preserve 
the probability density And it is straightforward to check the the wider Shapo- 

valov and Sukhomlin equivalence group breaks this rule, because it, generally speaking, 
does not preserve By this very reason, we restrict our considerations to the gauge 

transformations only. 

Definition 4 Two solutions of the Pauli equation with separated variables are called 
equivalent if they can be transformed one into another by group transformations from the 
Lie transformation group admitted by Pauli equation Moreover, solutions of the Pauli 
equation with separated variables having equivalent (in the sense of equivalence relation 
(^-(flZ^ spectral parameters X are equivalent. 



Definition 5 Two coordinate systems t, ui, u>2, uj 3 andt' , uj' 2 , uj' 3 are called equivalent 
if they give equivalent solutions with separated variables. In particular, two coordinate 
systems are equivalent if the corresponding Ansatze (fjj are transformed one into another 
by reversible transformations of the form 

t->t' = f (t), u a -> J a = f a (u a ), a = 1,2,3, (13) 
Q^Q' = QlomMkMhM, (14) 

where fo, ■ ■ ■ , f 3 are some smooth functions and Iq, . . . , l 3 are some smooth 2 x 2 -matrix 
functions of the indicated variables. 

Indeed, transformations (|13D and ([i~4]) preserve the form of Ansatze (||). So after com- 



pleting the procedure of separation of variables in these coordinate systems we obtain the 
same solutions with separated variables. 

These equivalence relations reflect the freedom in choice of the functions Q, u>i, u 2 , uj 3 
and separation constants Ai,A2,A3 preserving the form of the conditions (0)-®. They 
split the set of all possible vector-potentials, providing separability of equation (0), and 
sets of solutions with separated variables and corresponding coordinate systems into equiv- 
alence classes. In a sequel, when presenting the corresponding lists we will give only one 
representative for each equivalence class. 
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2 Classification of separable Pauli equations ([t]) 

In this section we obtain an exhaustive classification of the Pauli equations solvable within 
the framework of the approach described in the Introduction. Furthermore, we describe 
curvilinear coordinate systems enabling separation of variables in ([!]). 
Using the equalities (0) and (J6|) we get 

+ PuaXa] = 0. 

Splitting the expression with respect to A a yields 

[P fia ,Puf>] + [Ptf,Pu a ] = 0, (15) 

where /i, v,a,(3 = 0, 1, 2, 3 and henceforth summation over repeated Greek indices is not 
used. Choosing a = (3 we have 

LP P 1 = 

Taking into account this equality and the fact that any Hermitian (2 x 2)-matrix can 
be represented as a linear combination of the unit and Pauli matrices fl2|), we get the 
following form of P^ a : 

Pfia Pfj,a 

(o; M )J + G lxa {u fJ )s a a, (16) 

where F^, G m are some smooth scalar functions of the indicated variables, Uq = t and 
s a is a constant three-component vector. Substitution of expression ( |T6"D into (p~5|) yields 

{G^aGyp — G fj,pG V0 ^[s a a ', s^c] = 0. 

From this equality we conclude that there are two distinct cases: either s a ~ Sp or 
Gfj, a ~ G^/j. In view of this fact we get the two possible forms for the equations (0): 

z0 o = - (F 00 (t) + F ob (t)X b + (Goo(i) + G 06 (t)A fe )s<j) <^ , (17) 
<p a = (F a0 (uj a ) + F ab (u a )\ b + (G a0 (u a ) + G a b(v a )\b)sa) fa 

and 

i<Po = - (Foo(t) + F 6(t)A b + G (t)(s + s b X b )a) </? Q , (18) 
V^a = (F a0 (uj a ) + F ab (u a )X b + G a (u; a )(so + s b \ b )a) (p a , 

with a = 1, 2, 3. 

Definition 1 is quite algorithmic in the sense that it contains a regular algorithm of 
variable separation in Pauli equation ([I]). Formulas (|5|), ([XT|) — (p~8|) form the input data of 
the method. The principal steps of the procedure of variable separation in Pauli equation 
I 1 j £1X6 clS follows. 

1. We insert the Ansatz (|5|) into the Pauli equation and express the derivatives ifio, tp\, 
(f> 2 , <p 3 in terms of functions (p , <fi, <p 2 , f3, using equations flT7|)-(|r8"D. 
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2. We regard (fo, <pi, (p 2} f3, Xi, X 2 , X 3 as new independent variables y±, . . ., yj. As the 
functions Q, u>i, u 2 , u>s, A , Ai, A 2 , A 3 are independent on the variables yi, ■ ■ ■ ,2/7, 
we can demand that the obtained equality is transformed into identity under arbi- 
trary ?/!,..., y 7 . In other words, we should split the equality with respect to these 
variables under condition of commutativity @. After splitting we get an overdeter- 
mined system of nonlinear partial differential equations for unknown functions Q, 
ui, uj 2 , uj 3 , A , Ax, A 2 , A 3 . 

3. After solving the above system we get an exhaustive description of vector-potentials 
A(t,x) providing separability of the Pauli equation and corresponding coordinate 
systems. 

Having performed the first two steps of the above algorithm we obtain the system of 
nonlinear matrix PDEs: 

(i) |r|r = 0, 6/c, 6,0 = 1,2,3; 

E^M^^ = ^(t), 6 = 1,2,3; 

a=l c c 

For case of reduced equations given by ([17]) 
G wM dx a dx a = G 0lt (t), /i = 0,l,2,3; 

a=l c c 



III 



\111a) 



(iiib) For case of reduced equations given by (|n 

dug du c 

dXr dXr 



a=l 



Ur) 2(|^- ieQA b ) ^ + Q(t^ + Auj a \ = 0. a = 1.2. 3: 



dxb I dxh V dt 

3 



t \ n^sT t? t \ 9uJ aduj a .dQ dQ 

a=l 

+ (-F 00 (t) - - eA - e 2 A b A b + eaH^j Q = 0. 

Thus the problem of variable separation in the Pauli equation reduces to integrating of 
a system of nonlinear PDEs for eight unknown functions A , Ai, A 2 , A 3 ,Q,lji,uj 2 ,uj 3 of 
four variables t, x. What is more, some coefficients are arbitrary matrix functions which 
should be determined in the process of integrating of the system of PDEs (i) — (v). We 
succeeded in constructing the general solution of the latter which yields, in particular, all 
possible vector-potentials A(t, x) = (A (t, x), . . . , A 3 (t, x)) such that Pauli equation ([[]) is 
solvable by the method of separation of variables. 

In view of (ED we can always choose from each set of the equations (ii) — (iiia) and 
(ii) — (iiib) three such equations that the matrix of coefficients of 0J aXc ^>ax c (a=l,2,3) is 
non-singular. It is called the Stackel matrix [ITS]. The system consisting of these three 
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equations and of the equations (z) was integrated in 0. Its general solution u = u(t,x) 
is given implicitly within the equivalence relation by the following formulas: 

x = 0(t)£(t)(z(co)+v(t)), (19) 

Here 0(t) is a time-dependent 3x3 orthogonal matrix with Euler angles a(t), /?(£), 7(t): 

(cos a cos (3 — sin a sin (3 cos 7 — cos a sin f3 — sin a cos (3 cos 7 sin a sin 7 
sin a cos (3 + cos a sin (3 cos 7 — sin a sin /3 + cos a cos /3 cos 7 — cos a sin 7 
sin (3 sin 7 cos (3 sin 7 cos 7 

(20) 

v(t) stands for the vector-column whose entries t>i(£), v 2 (t), v 3 (t) are arbitrary smooth 
functions oft; z — z(u) is given by one of the eleven formulas 

1. Cartesian coordinate system, 

Z\ = LO\, Z2 = LO2, Z 3 = LO3, 

loi,uj 2 ,uj 3 e R. 

2. Cylindrical coordinate system, 

Z\ = e Ul cosu 2 , z 2 = e^ 1 sina; 2 , z 3 = U3, 
< uj 2 < 2n, u)i, u 3 G R. 

3. Parabolic cylindrical coordinate system, 
zi = (uj\ - u%)/2, z 2 = U1U2, z 3 = u 3 , 

Ld 1 > 0, Cg> 2 , UJ 3 G R. 

4. Elliptic cylindrical coordinate system, 

z 1 = a cosh c^i cos cl>2, z 2 = a sinh sin u 2 , z 3 = u 3 , 
lux > 0, —it < uj 2 < 7T, CU3 G R, a > 0. 

5. Spherical coordinate system, 
z x = cuf 1 sechco>2 cosco> 3 , 

z 2 = to i 1 sech o>2 sin co> 3 , 
z 3 = tuf 1 tanhcj2, 

> 0, w 2 e R, < u 3 < 2n. 

6. Prolate spheroidal coordinate system, 
z-i = a csch u>i sech u 2 cos 1^3, a > 0, 
z 2 = cl csch u;i sech ui 2 sin a; 3 , 

Z3 = acothcji tanhco>2, (21) 
uji > 0, w 2 e R, < u 3 < 2n. 

7. Oblate spheroidal coordinate system, 
Z\ — a esc uj\ sech u 2 cos u 3 , a > 0, 
z 2 = a esc uji sech a; 2 sin a; 3 , 

z 3 = a cot c^i tanhu^, 

< LOi < 7r/2, w 2 e R, < co> 3 < 2n. 
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8. Parabolic coordinate system, 
Zx = e" 1 ^ 2 coswa, z 2 = e^ 1+W2 sinw 3> 
z 3 = (e 2 ^-e 2 ^)/2, 

Wi,W2 G R, < uo 3 < 2tt. 

9. Paraboloidal coordinate system, 
Z\ = 2a cosh uJi cos uj 2 sinh u 3 , a > 0, 
z 2 = 2a sinh iO\ sin a; 2 cosh cu 3 , 
z 3 = a(cosh 2io\ + cos 2uj2 — cosh 2u 3 )/2, 

CJi, G Ft, < Co> 2 < 7T. 

10. Ellipsoidal coordinate system, 

Zi — a — f— — — dn(cj 2 , kf) sn(u 3 , k), a > 0, 
sn(u;i, k) 

da(ui, k) , 

z 2 = a — r- cn(w 2 , k ) cn(w 3 , k), 

sn.(ui, k) 

cnK, k) 

z 3 = a — — sn(u; 2 , k ) dn(^ 3 , k), 

sn(co>i, k) 

< u) x < K, —K' < uj 2 < K', 0<u 3 < AK. 

11. Conical coordinate system, 
Z\ = to , f 1 dn(cj2, k')sn(u! 3 , k), 
z 2 = w ] " 1 cn(co'2, k') cii(uj 3 , k), 
z 3 = LU^ 1 sn(cu 2 , k') dn(u; 3 , k), 

wi > 0, —K' <u 2 <K', 0<u 3 < AK; 

and C(t) is a 3 x 3 diagonal matrix 

h(t) 

C{t) = I l 2 (t) I , (22) 


where h(t) , l 2 (t) , l 3 (t) are arbitrary non-zero smooth functions that satisfy the following 
conditions 




• l\(t) = l 2 (t) for the partially split coordinate systems (cases 2-4 from (|2T|)). 

• li(t) = l 2 {t) = l 3 (t) for non-split coordinate systems (cases 5-11 from fl2l"|)). 

Here we use the usual notations for the trigonometric, hyperbolic and Jacobi elliptic 
functions, number k (0 < k < 1) being the modulus of the latter and k' = (1 — k 2 ) 1 ! 2 . 

From a geometric point of view the right-hand side of formula (|TP| ) is a result of 
application to vector z(u) of the following time-dependent transformations performed 
one after another: 

1. translations z ^ z 1 = z + v(t), 
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2. dilatations z — > z 1 = C(t)z, 

3. three-dimensional rotations z — > z* = 0(t)z with Euler angles a(t), Pit), j(t). 



Together with the rotations the following vector Q(t) = (Qi,Q 2 ^a) is considered [|T^, 
§35] 

Qi(t) = j(t) cosa(t) + $(t) sina(t) sin7(t), 

Q 2 (t) = 7(f) sin a(t) -$(t) cos a(t) sin 7(f), (23) 

fi 3 (t) = d(t) + /3(t) cos 7$, 

that is directed along momentary axis of rotation and called angular velocity vector. 

Note that we have chosen the coordinate systems u)i,u)2, 0)3 by means of the equivalence 
relation (|13|) in such a way that the relations 

Au a = 0, a = 1,2, 3 (24) 

hold for all the cases 1-11 in (|2T|) . 

After integration of system (i)—(iii) it is not difficult to integrate the remaining equa- 
tions (iv) and (v) from the system under study, since they can be regarded as algebraic 
equations for the functions A a (t,x), (a = 1,2,3) and A (t,x), correspondingly. 

Multiplying equation (iv) from the right by Q _1 we obtain for each component of matri- 
dQ 

ces — — Q 1 , b = 1, 2, 3 the systems of three linear algebraic equations. The determinants 

OXb 

of the systems do not vanish according to So, they have the unique solution 

^-Q- x = f b (t,x)I, 6=1,2,3, (25) 

OXb 

where fb(t, x) are scalar smooth functions and I is unit 2 x 2-matrix. From the compati- 
bility conditions 

9fa dfb , 100 

7T~ = x~ ' a ,6=l,2,3 
ox b OXa 

of the above system of PDEs we obtain that there exists such function g(t,x) that the 
equalities f a = dg/dx a , a = 1, 2, 3 hold. So (E51) takes the form 



^— = ^— Q, 6 = 1,2,3. 
ax;, ax;, 

The general solution of this system of matrix PDEs is 

Q = U(t) expg(t,x), (26) 

where U(t) is arbitrary 2 x 2-matrix function of t. 

Let us represent the complex-valued function g(t,x) in (|26|) as g = Si + iS, where 
Si, S are real- valued functions. Now, if we take into account that the components of 
the vector potential A(t,x) and functions 0)1,0)2,0)3 are real-valued functions, then after 
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inserting ( ptj| ) into (iv) with the use of fl2~3| ) we can split the obtained equations into real 
and imaginary parts: 

dS 1 du a 

^—^— = 0, a = 1,2, 3; 27 

OXb OXb 

( dS , \ du a doj a , , 

2 {Er b - eAb )ai + ~9f = °- a = 1 ' 2 ' 3 ' (28) 



Taking into account the equality @, we obtain from (|27|) the equalities dSi/dxb = 0, b = 
1, 2, 3. It gives that Si = Si(t). 

Let us denote eA = eA — VS. Then the system (p8|) takes the form of three linear 
algebraic equations for functions Al, A2, A3: 

du a duj a 

—- = 2e——A b , a =1,2, 3. 
at OXb 

The determinant of this system does not vanish due to @. Consequently, it has a unique 
solution. Making in this solution the hodographic transformation 

t = t, x a = u a (t,u 1 ,u 2 ,u)3), a = 1,2, 3, (29) 
we get the following expressions for A\, A2, A3: 

1 du(t,0) 



A 



2e dt 



After substitution into this formula expression for u(t,u) ([19|) , we return to variables 
t,Xi,X2,X3 and thus obtain the following system: 

2(-eA(t, x) + VS) = M (t)x + 0(t)C(t)il. (30) 

Here we use the designation 

M(t) = 6{t)0-\t) + 0(t)t(t)£r l {t)0-\t), (31) 



where 0(t), C(t) are variable 3x3 matrices defined by formulas (p0| ) and (P2|) , corre- 
spondingly, v = (vi(t), v<2(t), f 3 (t)) T . Note that 00~ x is antisymmetric and O tCr x O~ x 
is symmetric part of matrix A4. 

The direct calculation shows that equations (v) and (BO) are invariant under gauge 
transformations (|TTD. Thus the function S is transformed by the rule 

S ^ S' = S + ef, (32) 



which follows from (pa). In other words, if the transformations (11), (p2) in equations (v 



and (p0|) are performed one after another, we obtain the initial equations where functions 
A, A , S should be replaced with functions A', A' Q , S'. So, if the Pauli equation ([[]) with 
potential A, Aq admits separation of variables in some coordinate system, then the Pauli 
equation with potential A', A' admits separation of variables in the same coordinate sys- 
tem (the multiplier Q ( ^q) is changed only). Therefore, it is worthwhile to fix some gauge 
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and to work only with representatives of the equivalence classes of potentials A(t, x) (in 
the sense of equivalence relation (11)). 

We choose the gauge in a way that the equality 

2VS = 0(t)t(t)C- l {t)0-\t)$ + 0(t)C{t)'v. (33) 

holds. After integration of this system of PDEs we obtain the expression for S: 

s = \ E \r< + awxj > (34) 

where we use the notations 

x' = Q~ x x. (35) 

Next, we obtain from equation (|30|) the explicit form for space- like components of 
vector-potential of electromagnetic field 

A(t,x) = -—00- 1 x, (36) 

where the explicit form of matrix OO^ 1 is given by the formula 

(0 — (a + $ cos 7) 7 sin a — (3 cos a sin 7 

d + /5cos7 — (7 cos a + (3 sin a sin 7) (c]7) 

— (7 sin a — /3 cos a sin 7) 7 cos a + j3 sin a sin 7 

where a, /3, 7 are arbitrary functions of t. 

Thus formula ( ^6|) means that the space-like components of electromagnetic field A(t, x) 
are linear with respect of spatial variables. So the magnetic field H = rot A should be 
homogeneous, i.e., independent of spatial variables x. From formulas (0), (|37|) we can 
obtain its explicit form 

eHi = —j(t) cos a(t) — $(t) sina(t) sin7(t), 
e H 2 = -7(t) sina(t) + f3(t) cosa(t) sm^y(t), (38) 
eif3 = -«(*) - Pit) cos7(t). 

Now the space-like components of the electromagnetic field take the final form 

! / -ff 8 (*) ^ 2 (t) \ x 
A(*,f) = - F 3 (t) -H^t) \x = -H(t)xx, (39) 
2 V -H 2 (t) H x {t) / 2 

where symbol x denotes cross product. 

Within the equivalence relation (|HD we can always choose the function ti (t) to be a 
solution of matrix ODE 

iU = {-eaH(t))U (40) 

with the initial conditions W(0) = /. Due to the theorem of existence and the uniqueness 
of the solution of the Cauchy problem for the system of ODEs there is unique solution 
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U{t) of system (HO) for each fixed configuration of magnetic field H(t). Moreover, matrix 
U{t) is a unitary one. Indeed, taking into account we have the equality 

j t (U*U) =U*{ieaH)U +U*(-ieaH)U = 0, 

i.e., U*U = const. The initial conditions give U*U = I. 

Thus we can consider the following change of variables in the Pauli equation ([!]) 

i) = u{t)i). (41) 

Due to the unitarity of the matrix U the quantity which is regarded in quantum 
mechanics as the probability density, is not changed. So the change of variables ([IT]) is the 
correct one. As a result, the term eaH in Pauli equation (HD vanishes, and we obtained a 
system of two Schrodinger equations for the function ip. 
Thus we proved the following assertion. 

Lemma 1 A necessary condition for the Pauli equation (fjj to be separable (in the sense 
of definition is that it has to be equivalent (in the sense of equivalence relation ^D) 
to a system of two uncoupled Schrodinger equations. 

Let us substitute the equality ( |26l) into equation (t>), taking into account equations 



(flOl) and S\ = Si(t). Splitting the equation obtained into real and imaginary parts (note 
that all functions F 00 , F a0 , a = 1, 2, 3 are real-valued ones), we obtain the equalities 

E„ , .du a duj a dS OS OS . OS , . . , . . . . 

i^ooK ) ^ ^ -^r- T^TT" + 2eA fe — - F 00 (t) - eA - e 2 A b A b = 42 
oxh ox b at ax b ax b ax b 

a=l 

Si + AS - e^P- = 0. (43) 
axh 



Inserting into equation (J42[) expressions for S (gj) and Ai,A 2 ,A 3 fl39|) , we obtain the 
explicit form of Aq. 

3 

eA (t, x) = J2 F -oM^^ ~ F 00 (t) - e 2 A b A b - I P. (44) 

a=l ^ ^ 



Here A b A b follows from (|39]), @ 

AA b A b = (H 2 x 3 - H 3 x 2 ) 2 + {H 3Xl - H lX3 ) 2 + (H 2 x 1 - H lX2 ) 2 , (45) 
where Hi,H 2 , H 3 are components of magnetic field ([38D; function P has the form 



3 



a=l 



la 



where x[,x 2 ,x' 3 are given by formula (j35l) and l a = l a (t), v a = v a (t), a = 1,2,3, are 
arbitrary smooth functions, which define new coordinate system fll9|) . 
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Let us emphasize that the expression for Aq includes arbitrary functions Fio(ux), 
-^20(^2), F 3 o(uj 3 ) , F 00 (t) , where functions u a = u a (t,x), a = 1,2,3, belong to one of 11 
classes, whose representatives are given implicitly by the formulas fli~9l) -( p^) . 

Below we give explicit forms of the eikonals R~ 2 = -rr—^r~ f° r eacn class of ui a (see, 

OXb OXb 

also M): 



1. 


Rf 


= K 2 , 


i = 1,2,3; 


2. 


Ri 2 


= R 2 2 -- 


= K 2 e- 2 -\ Rf = hf; 


3. 


RI 2 


= R 2 2 - 


= hf(co 2 1 +u J 2 2 )-\ Rf = hf; 


4. 


Rf 


= R 2 2 -- 


= h^ 2 a~ 2 (cosh 2 lu 1 — cos 2 cl> 2 ) -1 , -R7 2 = h^ 2 ; 


5. 


Ri 2 


= h± 2 uj 


f, R 2 2 = R 3 2 = K[ 2 <J\ cosh 2 ^; 


6. 


K 


= h\ 2 a 


" 2 sinh 2 ^(sinh -2 uj\ + cosh -2 u^)^ 1 ? 




R2 


= h\ 2 a 


" 2 cosh 2 UJ2 (sinh -2 uo\ + cosh -2 cl^) -1 , 




R3 2 


= h\ 2 a 


" 2 sinh 2 uji cosh 2 cu 2 ; 


7. 


K 


= h\ 2 a 


~ 2 sin 2 uji (sin -2 uj\ — cosh -2 u^) -1 , 




R 2 


= K{ 2 a 


" 2 cosh 2 co> 2 (sin" 2 uji — cosh -2 c<-> 2 ) -1 , 




R3 2 


= h\ 2 a 


" 2 sin 2 u>i cosh 2 tu 2 ; 


8. 


R \ 


= hfe 


~ 2uJl (e 2 ^ +e 2 " 2 )-\ 




R 2 


= hfe 


' 2 ^(e 2uJl + e 2uJ2 )-\ Rf = h^ 2 e~ 2(uJl+UJ2) ; 


9. 


R \ 


= h\ 2 a 


~ 2 (cosh 2u>i — cos 2co>2) -1 (cosh 2uj\ + cosh2u;3) - 




R 2 


= h^ 2 a 


~ 2 (cosh2co>i — cos 2kj 2 )~ 1 (cos 2cc>2 + cosh2co> 3 ) -1 . 




Rf 


= h\ 2 a 


~ 2 (cosh2c(Ji + cosh 2cu 3 )~ 1 (cos 2t02 + cosh2c<j3) - 


10. 


Ri 2 


= h\ 2 a 


-2 ( dn 2 (co>i, k) 2 2 X" 1 /dn 2 (wi, A;) 

U 2 K*) kcn ^ k) ) U a K*) 



(47) 



+ fc 2 cn 2 (c<j3, k) 



R 2 2 = K 2 a- 2 ( ^^1 - ^cn 2 ^, fc')) (^cn 2 ^, if) + fcWfa, 

i?- 2 = hfa- 2 ( ^1^1 + fc 2 cn 2 (^ 3 , k)) 1 (£/ 2 cn 2 (u, 2) k') + k 2 cn 2 (u 3 , k))' 1 ; 
\sn 2 (uji,k) J 



11. i?7/ 2 = hfuf, R 2 2 = Rf = h^ 2 uo\ (A; /2 cn 2 (w 2 , k!) + A; 2 cn 2 (w 3 , k)) 1 . 



At last, let us find the multiplier Q. Substituting the formulas ( p4|) and fl39|) into 
equation fl43| ) gives 

1 3 /' 

a=l 

whence it follows that 

1 3 

S! = --X>1„. (48) 
0=1 

Taking into account expression for 5 (|34]) , we obtain from formula (|26|) the explicit form 
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of Q 

Q = U(t)-^=exp g I (^-x'l + 2l a v a x'\ , (49) 

where U(t) is given by the equation (fiO|), and a^, x 2 , x' 3 are given by formula (|35|) . 
Thus we have proved the main result of the article: 

Theorem 1 Pauli equation ([][) admits separation of variables (in the sense of definition 
[![) if and only if it is gauge equivalent to Pauli equation where 

- the magnetic field H = rot A is independent of the spatial variables, 

- the space-like components A\, A 2 , A 3 of the vector-potential of the electromagnetic 
field are given by (fM 



- the time-like component Aq is given by formulas (f^j-ffify. 

Comparing the components of magnetic field fl38|) with components of angular velocity 
vector ( ^3|) of rotation of coordinate system fll9l), we obtain the equality eH = —Q. So, 
we prove the following assertion: 

Corollary 1 Let Pauli equation (fl[) admit separation of variables in some non- stationary 
coordinate system t, u a = u> a (t,x), a = 1,2,3, where functions ui(t,x), co> 2 (t,x), u 3 (t,x) 
are given implicitly by formulas ^T^j-j^j). Then angular velocity vector ft2$j ) of rotation 



of this coordinate system equals —eH , where H = rot A is magnetic field. 

It follows from the corollary that a necessary condition for the Pauli equation (0) with 
non-zero magnetic field H to be separable (in the sense of our definition 1) is that the 
angular velocity vector ( p3|) of rotation of the separation coordinate system (|i9|) -(p2"l) has 
to be non-zero. 

Summing up we conclude that coordinate systems and vector-potentials of the elec- 
tromagnetic field A(t,x) = (A (t , x) , A(t , x)) providing separability of the corresponding 
Pauli equations coincide with those for the Schrodinger equations. Namely, we prove 
that the magnetic field H = rot A has to be independent of the spatial variables. Next, 
we have eleven classes of potentials A (t,x), corresponding to eleven classes of coordi- 
nate systems t,u a = u a (t,x), a = 1,2,3, where the functions ui(t, x), u>2(t, x), uj 3 (t, x) 
are given implicitly by formulas (|i~9D-(p2j). Pauli equation ([I]) for each class of the func- 
tions A (t, x), A(t, x) defined by (|39|) , ([44]) and ( ^7|) under arbitrary F o(t), -Fao(^a) and 
fixed arbitrary functions a(t), Pit), j(t), v a (t), l a (t), a = 1,2,3, separates in exactly one 
coordinate system. 



The solutions with separated variables are of the form (|5|), where Q is given by (f49|) . 
The separation equations read as (p~7|) or (0), where the functions F^ , /i = 0, 1, 2, 3, are 
arbitrary smooth functions defining the form of the time-like component of the vector- 
potential A(t, x) (see, (|44|)). The explicit forms of other coefficients F^ a ,G ^ U ,G M of re- 
duced equations can be obtained by splitting relations (ii) and (Hi) with respect to inde- 
pendent variables ui,u 2 ,uj 3 ,t for each class of the functions z — z(lj) given in (^TJ). Let 
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us denote 



S 



( T X 

Sn 
S21 
V S 3 i 



T 2 
S12 
S22 

5*32 



^3 \ 

513 
5*23 
S33 / 



(50) 



where the functions S a b{uj a ) (a,b = 1,2,3) are given below as entries of 3 x 3 Stackel 
matrices, whose structure is determined by the choice of the functions z — z(uj): 



F A 



F & = 



F» 




a 2 cosh 2 ui 
—a 2 cos 2 u>2 




1 
-1 
1 



3 2wi 







F, 



lift 




-UJ-, 



F, 








COsh Cg> 2 





a 2 sinh 4 oj\ — sinh 2 uj\ — 1 
a 2 cosh - UJ2 cosh - ui2 — 1 
] 







F 



10 




F 9 



F 7 = 



a 2 cosh 2 2uji 
—a 2 cos 2 2u>2 
a 2 cosh 2 2uz 



u sin Ui — sin 2 oj\ 1 

a 2 cosh -4 u;2 cosh -2 co> 2 — 1 




2 • —4 

a sm 

-n 2 



—a cosh 2oj\ 
a cos 2u>2 
a cosh 2u;3 




dn (u>x, k) 
sn 4 (c<Ji, k) 
-a 2 k' 4 cn 4 (u 2 ,k') 
\ a 2 k 4 cn 4 (u;3, k) 



dn 2 (c<Ji, k) 

sn 2 (a;i, k) 
k' 2 cn 2 (u 2 ,k') 
k 2 cn 2 (cj 3 , k) 



\ 



1 / 



F 
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k 2 cn 2 (u) 3 , k) 





-1 
1 



The functions Ti(t),T2(t),T^(t) are expressed in terms of the functions hi(t), h2(t), hs(t): 

(52) 



1. 

2-4. 
5- 11. 



Ti 
Ti 
T-i 



K 2 
,-2 



% = 1,2,3; 
h^, T 2 = 0, T 3 
hf, T 2 = T 3 = 0. 



Let and M be 3 x 3 constant matrices. Now, if the reduced equations are given by 
(0), then 

rp II j? ||3 3 ||^ m3 3 

~~ 11^^11^=0 0=1' ^ — 11^^11^=00=1 

are block (6 x 8)-matrices, where F^ a and G^ a are 2 x 2-matrices that are equal to products 
of the corresponding entries of the matrices SK and SM by the unit (in the case of the 
matrix F) or sa (in the case of the matrix G) matrices. Accordingly, equation (§) takes 
the form 

rank (F + G) = 6. (53) 

If rank K = 3, then we can always rearrange Ai, A2, A3 with the use of the equivalence 
relation (H) in order to get K = I. Analogously, without loss of generality we may put 
M = I, provided rankM = 3 and s 2 7^ 0. 
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1 1 1 1 3 — * — * 

If rank M = 0, then the column ||G>o|| M= o has necessarily the form Sg, where g is a 
constant three-component column. If rank M^O, then we can always kill this column 
by a proper rearranging of Ai, A2, A3 with the use of the equivalence relation (||). 

Next, if the reduced equations are given by (|IBD , then the matrix F is defined in the 
same way as in the previous case. Furthermore, 

G = WGpSaffW^^ 

is a block (6 x 8)-matrix, where are the three-component columns Sg^ (g^ is a constant 
three-component column). In addition, in this case identity Q55p holds, so that we can 
put K — I, when rank K = 3. If s a , (1 = 1,2,3) are three linear independent vectors, 
then we can always put Sq = 0. 

We will finish this section with the following remark. It follows from theorem 1 that 
a choice of magnetic fields H allowing for variable separation in the corresponding Pauli 
equation is very restricted. Namely, the magnetic field should be independent of spatial 
variables xi,x 2 ,x 3 in order to provide the separability of Pauli equation (|l|) into three 
second-order matrix ordinary differential equations of the form (|7|). However, if we allow 
for separation equations to be of a lower order, then additional possibilities for variable 
separation in the Pauli equation arise. As an example, we give the vector potential 



A(t, x) = \^Ao \ Jxl + xlj , 0, 0, A 3 \ Jxl + x 

where A , A 3 are arbitrary smooth functions. The Pauli equation (jH) with this vector-po- 
tential separates in the cylindrical coordinate system 



t, u>i = In x\ + x\j , uj 2 = arctan(a;i/a;2), ^3 = £3 

into two first-order and one second-order matrix ordinary differential equations. The cor- 
responding magnetic field H = rot A is evidently x-dependent. In this respect, let us also 
mention the recent paper by Benenti with co-authors 0, where the problem of separation 
of variables in the stationary Hamilton- Jacobi equation with vector-potential has been 
studied. They have presented a number of vector-potentials, for which the Hamilton- 
Jacobi equation is separable, and the corresponding magnetic fields are inhomogeneous 
ones. These potentials allow for separation of variables in the stationary Schrodinger 



and Pauli equations with vector-potentials as well (see, e.g., |15j concerning the relation- 
ship between the separation of variables in the Schrodinger and Hamilton- Jacobi equa- 
tions). These facts imply an importance of application of our approach to classify the 
non-stationary Pauli equations of the form (|I|), which admit separation of variables into 
first- and second-order matrix ordinary differential equations. We remind that here we 
give the classification results for the case, when all the reduced equations are second-order 
ones. We intend to address this problem in one of our future publications. 
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3 Algorithm of separation of variables in the Pauli 
equation with fixed potential 

Theorem [I] gives the solution of the problem of classification of the Pauli equations ([I]) 
with variable coefficients that are separable (in the sense of definition 1) at least in one 
coordinate system. 

Let us consider the problem of classification of coordinate systems that allow for sep- 
aration of variables (in the sense of definition 1) in the Pauli equation (jl]) with fixed 
vector-potential A , A. 

Let some fixed vector-potential A(t,x), Ao(t,x) be given. The scheme of finding all 
coordinate systems providing separation of variables is as follows: 

1. With help of gauge transformations (|ll]) we reduce the space-like components of 
vector-potential A(t,x) to the form (|39|) . If it is impossible, then Pauli equation ([!]) 
with this vector-potential is not solvable by the method of separation of variables 
in the framework of our approach. 



2. We solve the system of ODE ( |3"%D for given magnetic field H(t) and obtain the 
explicit form of functions a(t), j3(t), 7(t). 

3. For each of 11 classes of coordinate systems t, u a = u a (t,x), a = 1,2,3, which are 
given by formulas (|I~9|)-(p2j), taking into account restrictions obtained on the first 



step of the algorithm, we find the explicit form of 

(a) the time-like component A of the vector-potential in terms of ou; 

(b) function P, substituting in ( ^6[ ) the expression for x' in terms of a; (see formulas 



and ([ID): 

x' = C(t) (z(uj) + v(t)) ; (54) 
quantity e 2 A b A b by the formula 

Ae 2 A b A b = (n 2 x' 3 - n s x 2 ) 2 + (n 3 x[ - n^) 2 + (n 2 x[ - n x x' 2 ) 2 , 

where x[,x 2 ,x' 3 are given by the formula ([54"D, and functions rii,n2,n 3 are as 
follows 

rii = j(t) cos/3(t) + a(t) sm(3(t) sm^ft), 

n 2 = sin/3(t) + a{t) cos f3{t) shi7(t), (55) 

n 3 = $(t) + a(t) cos7(t); 

dx b dx b 
for given class of coordinates. 



(d) eikonals ^li^fl = R 2 a = 1,2,3, which are determined from the list (|47|) 

ox b ox b 



4. We substitute the equalities obtained into equation ([P]) and obtain 11 equations for 
each of 11 classes of coordinate systems t, u>i, u 2 , u 3 . For each of these equalities we 
find all possible functions F a0 (u; a ), a = 1,2,3, F 00 (t) that reduce it to the identity 
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by the independent variables t, ui, u 2 , uj 3 (i.e. we split this equality with respect 
to these variables). It gives, in its turn, the explicit form of the functions v a (t), 
l a (t), a = 1,2,3 and additional restriction on a(t), /3(t), j(t), giving the form of 
the coordinate system in question. All obtained coordinates for which the functions 
F a o(uj a ), a = 1,2,3, F 00 (t) exist are only coordinate systems providing separability 
of Pauli equations in the sense of definition [TJ. 

Example. As illustration of this algorithm consider the problem of separation of vari- 
ables in Pauli equation for a particle interacting with a constant magnetic field. With- 
out loss of generality we can always choose it as directed along axes OZ: eH = (0, 0, c) T , 
where c is a non-zero real constant. The vector-potential of electro-magnetic field has the 
form 

2eA = c Of, eA = p- {x\ + x\ - 2xf) , (56) 

\ / \x\ 12 

where q is a non-zero real constant and \x\ = \/ xf + x\ + x\. 

A direct check shows this vector-potential satisfies the vacuum Maxwell equations 
without currents 

DA -^- f^+divA) =0, (57) 



dt V dt 
uA + grad + div A^j = 0, 

where □ = d 2 /dt 2 — A is d'Alembert operator. Therefore, it is a natural generalization 
of the standard Coulomb potential, which is obtained from (|56|) under c — > 0. 

Proposition 1 The set of inequivalent coordinate systems providing separability of the 
Pauli equation ([][) with vector potential of electromagnetic field ( [56[ j is exhausted by the 
following ones: 

x = 0(t)z, (58) 



where O is a time- dependent 3x3 orthogonal matrix t^Uj), with Euler angles 

a{t) = —ct, (3 = const, 7 = const, (59) 
and z is one of the following coordinate systems: 
1. spherical (formula 5 from list ft2l\)), 



2. prolate spheroidal II (formula 6 from (pH^, where one should replace Z3 with Z3 
a(coth oj\ tanha; 2 ± 1) ), 



3. conical (formula 11 from (\21\)). 

Proof. The space-like component A(t,x) of the given vector-potential (|56|) is already 
reduced to form (|39T). 
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The system of ODE ( |5SD for given magnetic field takes the form: 

7 cos a + (3 sin a sin 7 = 0, 7 sin a — /3 cos a sin 7 = 0, a + (3 cos 7 = — c. 

This implies the equivalent system 

7 = 0, /5sin7 = 0, d + /?cos7 = — c. 

Its general solution up to translation by t is given by formulas ( |55| ) (solution a ± (3 = —ct 
for case sin 7 = is included into (|59|) as a particular case after denoting a ± /? — > a) . 

The steps 3 and 4 of the above algorithm will be illustrated by the case of spherical 
coordinate system 5 from fl2"T|) (for other coordinate systems this procedure is an analogous 
one). For this case the equality (|44]) in terms u>i, u 2 , u 3 takes the form 

l~ 2 (FxoMut + (F 20 (uj2) + FsoMVi'cosh 2 ^) - F 00 (t) = 

3 

|ff + ( 2 (^ + 2^«K + /2?) a) , (60) 
0=1 

where Z = Zi = 1% = h, I 7^ (because of the spherical coordinate system is non-split one), 
and 

x[ = /(cuf 1 sechco> 2 cosco>3 + Vi(t)), 
x' 2 = l^u^ 1 sech.u 2 smu 3 + v 2 (t)), 
x' 3 = l^i 1 tanho^ + v 3 (t)). 

Next we perform on both parts of equality (BH) the following step by step operations: 

1. multiplying by lux, 

2. differentiation with respect to u>i, 

3. division by uj 2 , 

4. differentiation with respect to Mi, 

5. differentiation with respect to u 2 , 

6. multiplying by Ioj\\x' | 7 , 

7. twice multiplying by uj\. 
As result we get 

— 24g sech 3 a;2 (f 1 + v \ + v%)(vi cos 1^3 + v 2 sin 0^3 + t> 3 sinha^) x 
(— v 3 + (f 1 cos u; 3 + v 2 sin cj 3 ) sinh co> 2 ) = 0. 

The equality obtained is transformed into an identity with respect to independent variables 
lux, u>2, CJ3 if and only if the condition v± = v 2 = V3 = holds. Now the equality ( |60D takes 
the form 

r 2 (F 10 (^K + (F2o(^2) + F m (u 3 ))tu 2 1 cosh 2 u 2 ) -F 00 (t) = ^uj 1 + ij + lj) ( 61 ) 



Performing on both parts of equality (jBTf) the following step by step operations: 
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1. differentiation with respect to uji, 

2. multiplying by I 2 , 

3. differentiation with respect to t, 

4. multiplying by uf, 

5. differentiation with respect to oui, 

we get the equality 3qluf = 0. This implies / = const and with the help of dilatations we 
can put without loss of generality 1 = 1. Thus the coordinate system takes the form 
The equation (|Jl|) yields 

c 2 

Fio^i)^ + (F20O2) + F 30 (uj 3 ))uj 2 cosh 2 uj 2 - F 00 (t) = quj x + —uj x 2 . 

o 

We can split the equation obtained by the independent variables u)\, ui 2 , oj 3 . As a result 
we get 

c 2 

F10 = quj 1 3 + — 6 + hLU-L 4 - fc 2 ^i 2 , 


F 20 = fc 2 sech 2 u; 2 - k 3 , F 30 = k 3 , F 00 = k v 
The theorem is proved. > 



Separation of variables in the Pauli- Maxwell sys- 
tem 



The expressions (39), (Q)-(47) give the most general form of the vector-potential of 
the electromagnetic field, providing separability of the corresponding Pauli equations. 
But, because of generality of the results, these expressions are too cumbersome, and their 
physical interpretation is somewhat difficult. Therefore it would be interesting to know the 
form of these potentials under certain physical restrictions. The most natural restriction 
is that the vector-potential satisfies the vacuum Maxwell equations without currents (|57|) . 
In this section we describe all explicit forms of the vector-potentials A(t, x) that 

a) provide separability of Pauli equation, 

b) satisfy vacuum Maxwell equations without currents ( |o"7| ) and 

c) describe the non-zero magnetic field. 

Furthermore, we construct inequivalent coordinate systems enabling us to separate vari- 
ables in the corresponding Pauli equation. 

The similar problem with more strong restrictions was analyzed in [[?tj for a two-dimen- 
sional Schrodinger equation with vector-potential. Note that an analogous problem for 



the Dirac equation for an electron was analyzed in |2T . 
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Taking into account the form of A (^), the Maxwell equations ([)7|) take the form 

AA = 0, (62) 

and 

d 2 A d 2 A _ d 2 A 

— — — — — L3X2 + I2X3, — 7. — — h^i — hXi, 77^ — — —h^i + nX2- 

otoxi 010x2 010x2, 

From the compatibility conditions of the above system of PDEs we get 

h = i 2 = i 2 = 0, 



dtdx a 



0, a = 1,2, 3. (63) 



Inserting expression for potential Ao(t,x) ( |^ ) into ( |62"D with subsequent change of 
independent variables (|l^) yields (we use the relations Auji = 0, u)i Xa ujj Xa — 0, % ^ j, 
i,j = 1,2,3) 



£ i^)it7 2> ) ^ = jtr + e W + ^ 2 + H ^ (64) 



=1 * 



Rf= a 1 n * , z=l,2,3 (65) 




where the eikonals 



<9a; a <9z a : 
are given in the list (f4Tp. 

Thus we get eleven functional relations Pi, ... , Pn for each class of coordinate system 
(|19l), whose form is determined by the form of one of the eleven expressions z\(ui, u 2 , UJ3), 
z 2 (ui, u 2 , 003), zs(ui, u 2 , U3) from the list (f2ll). As t, Ui, uj 2 , 0J3 are functionally indepen- 
dent, we can split the above relations with respect to the variables t, Ui, u 2 , U3, thus get- 
ting ordinary differential equations for the functions F i0 (uji), k(t),i = 1, 2, 3. After solving 
them the formula ( [44]) yields the expressions for Aq in terms of variables t, oj\, u 2 , 003. 
Returning to variables t, xi,x 2 ,X3 (with the aid of (|19|) ) , we should split the expression 
obtained for A (t, x) with respect to t. Indeed, the general solution of the equation (|63|) 
is 

A (t,x) = f 1 (x) + f 2 (t). 

At the expense of the gauge invariance of the Pauli equation we may choose f2{t) = 
0. Thus the potential Aq should be a function of x only. This condition restricts the 
choice of A , thus giving ordinary differential equations for the functions kit) , Vi(t) , % = 
1,2,3. Solving them we obtain the explicit forms of the function -Fnn(£) and coordinate 
systems (0). After simplifying these coordinate systems with the aid of equivalence 
transformations we get a full description of the vector-potentials A(t, x) and coordinate 
systems, giving the solution of the problem under study. 

Omitting the details of the calculations (they are very cumbersome) we present below 
the results. Note, when presenting lists of the vector-potentials A(t, x) and coordinate 
systems we use invariance of the system of the Pauli and Maxwell equations with respect 
to the groups of rotations by spatial variables xi,X2,X3 and translations by all variables 
t,Xi,x 2 ,x 3 (see, e.g. 
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1. Case of non-stationary magnetic field: 

eH= (0,0, At + B), 

eA = — - {x\ + x\ - 2x\) + a 1 x 1 + a 2 x 2 + a 3 x 3 , 

where A, B, k, a%, a 2 , a 3 are arbitrary real constants. 
The coordinate system is 

x = CO(z + v). 

Here O is a time-dependent 3x3 orthogonal matrix 0(a, [3, 7), where 

a = —At 2 -Bt, (3 = 0, 7 = 0; 

z is Cartesian, cylindrical or elliptic cylindrical coordinate system (formulas 1, 2, 4 from 
(pl|); £ is the 3 x 3 diagonal matrix 

/ l(t) 
C=\ l(t) 

V i 3 {t) 

and v{t) is vector-column v{t) = (t>i, t>2, f 3 ) T where functions l(t),l 3 (t),vi(t),V2(t),v 3 (t) 
are solutions of the following system of ordinary differential equations: 

lv\ + 2lv\ + 4c-^j- — 2c u \ = — 2(ai cos a + a 2 sin a), 
r t 

v 2 1 

lii2 + 2/U2 + 4c— — 2ci2- = — 2(— ai sin a + a 2 cos a), 
I 6 I 



l 3 v 3 + 2l 3 v 3 + 4c 3 -3 - 2ci 3 — = -2a 3 . 

<o ta 



¥ " 2ci3 r 3 

Here c, c 3 , Cn, Ci 2 , C13 are arbitrary real constants. 
2. Cases of stationary magnetic field. 



Case 1: 



eH = (0, 0, k), k = const ^ 0; 

eA = ~y^( x i + x l ~ 2x 3) + a i x i + a 2^2 + a 3 ^3, 



where a = (ai,a2,a 3 ) is constant vector. 
The coordinate system is 

x = lO(z + v). 
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Here O is a time-dependent 3x3 orthogonal matrix 0(a,(3, r y), where a = —kt,(3 = 
const, 7 = const; z is one of coordinate system, given by formulas 1-11 from ([H]); function 
l(t) is solution of the equation 



2 3 I c 
k + 2~1 = ¥ 



given by one of the formulas: 



r = T l. f 2 = \/^±lsin [2sl\kt\ ( 



i- 



for coordinate system z given by the formulas 1, 2, 4, 5, 6, 7, 10, 11 from the list ( pTJ ) and 

c = 0, / = Ci sin | \f^kt 

for coordinate system z given by the formulas 1-11 from the list fl2T|) . Here Ci is an arbi- 
trary real constant. Vector v is a solution of the following system of ordinary differential 
equations: 

— 2c 

3/i; + 6ft> + — w = -dO~ l a. 
I 6 

Case 2: 

e# = (0, 0, k), k = const ^ 0; 

eAn = — i, + i, — 2xo), a = const ^ 0. 

The coordinate system is 

x = Oz. 

Here O is a time-dependent 3x3 orthogonal matrix 0(a, /3, 7), where a = —kt,/3 = 
const, 7 = const and is one of the following coordinate systems: 

1. spherical (formula 5 from (pl|)), 

2. prolate spheroidal II (formula 6 from (]2T|) , where one should replace 23 with 2:3 = 
a(cothu;i tanhu;2 ±1)), 

3. conical (formula 11 from (pl|)). 
Case 3: 

eH = (0, 0, k), k = const ^ 0; 

eX = -^(x? + xI-2xi) + ^ +O2 ?| + i|('|2ln^-l 
12 r r d r z \2r r — X3 

where r = \J x\ + x\ + rr§ and ai, 02, % are real constant numbers. 
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The coordinate system is 

x = lOz. 

Here O is a time-dependent 3x3 orthogonal matrix 0(a, j3, 7), where a = —kt, [3 = 7 = 0; 
z is the spherical coordinate system, given by formula 5 from (|2l|) and function l{t) is 
given by 

l 2 = ^C 2 ±^ sin U^toJ + Ci, or l = C 1 ^U^kt 

under condition a\ = and / = 1 under condition a± 7^ 0. Here C\ is an arbitrary real 
constant. 

Case 4: 

eH = (0, 0, k), k = const 7^ 0; 

. A; 2 , o 2 n 2\ a i a 2 / 1 , xt 1 . iA 

eA = — — [x\ + x 2 — 2xn) H H h a 3 — arctanh— arctanh — , 

^+ ^ \ j^~r <p~r np np / 

where xf = x 3 ± a and r ± = ^Jx\ + x\ + (x% ± a) 2 , and a,ai,a 2 ,a 3 are arbitrary real 
constants. The coordinate system is 

x = Oz. 

Here O is a time-dependent 3x3 orthogonal matrix 0(a, (3, 7), where a = —kt, (3 = 7 = 
and z is a prolate spheroidal coordinate system, given by formula 6 from (pi]) . 

Case 5: 

eH = (0, 0, k), k = const 7^ 0; 

eA = (x 2 + a; 2 , - 2x 2 ) + 2axa^ + 2a 2 -^ 2a 3 ( a^- arccot/i - arctanh-^- ) 

J-2 J JJi \ J JJi a JiJ 

where 



f = ^1 (a 2 - r 2 ) 2 + 4a 2 xl fi = \j ^2 ' r = V x i + x l + x §> 

and a, ai, a 2 , 03 are arbitrary real constants. The coordinate system is 

x = Oz. 

Here O is a time-dependent 3x3 orthogonal matrix 0{a, j3, 7), where a = —kt, [3 = 7 = 
and z is an oblate spheroidal coordinate system, given by formula 7 from fl2T|) . 
Note that expression for A Q can be rewritten in the form 

k 2 2 2 2\ a l + ^ fl 2 a l — 2<2 2 / 1 st 1 Xn 

eA = (x 1 + x 2 — 2x 3 ) -\ — 1 — h m 3 I — arctanh— — — arctanh- 



where xf = x 3 ± ia and f ± = ^x 2 + x\ + (x 3 ± ia) 2 . 
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Case 6: 



eH = (0, 0, k), k = const ^ 0; 

eA = — —ix 1 + x 2 - 2xn) H h a 2 x 3 H In , 

o r r r — X3 

where r = i/xf + x\ + x| and oi, 02, 03 are arbitrary real constants. 
The coordinate system is 

x = Oz. 

Here O is a time-dependent 3x3 orthogonal matrix 0(a, (3, 7), where a = —kt, (3 = 7 
and ifis a parabolic coordinate system, given by formula 8 from (21). 



Case 7: 

e# = (0,0, A;), k = const ^ 0; 

eA = - - (aci + X2 - 2x 2 3 ) + a m ( x i + ^2) + 0-3X3, 

where k, a, (13 are arbitrary real constants. 
The coordinate system is 

xi = e Ul cos(co>i — kt), x<i = e ui sin(c<Ji — kt), £3 = I3U3 + V3, 

where h,V3 are solutions of the system of ordinary differential equations 



5 _ 7T = ^ 3 + 2 ^ 3?)3 + 4C3 W _ 2c 137 L = -2^3- 
i 3 4 {3 i 3 {3 



Note that some of the potentials obtained have the clear physical meaning. For in- 
stance, cases 2 and 3 under condition k = a 2 = 03 = give the standard Coulomb 
potential. Case 4 under condition k = a 3 = gives the potential for a well-known two- 
center Kepler problem, i.e., the problem of finding wave functions of electron moving in 
the field of two fixed Coulomb centres with charges a%, a 2 and intercenter distance 2a (the 
model of ionized hydrogen molecule). Coulson and Joseph |2t| showed that the corre- 



sponding Schrodinger equation admits separation of variables in the prolate coordinate 
system only. We obtained this potential as a particular case of the more general potential. 



5 Concluding Remarks 

Theorem 1 provides the complete solution of the problem of classification of the Pauli 
equations ([!]), which are solvable within the framework of the method of separation of 
variables in the sense of our definition 1. According to these theorems the coordinate 
systems and the vector-potentials of the electromagnetic field A(t,x) = (A (t,x), A(t,x)) 
providing separability of the corresponding Pauli equations coincide with those for the 
Schrodinger equations with vector-potential. So the results obtained in the article are 
valid for the Schrodinger equation as well. 
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It is well-known that the possibility of variable separation in a system of PDEs is closely 
connected to its symmetry properties 0, fLOp . Namely, solutions with separated variables 



are common eigenfunctions of three matrix mutually commuting symmetry operators of 
the equation under study. For all the cases of variable separation in Pauli equation ([I]) 
such matrix second-order operators can be constructed in the explicit form, by analogy 
to what has been done in for the (l+2)-dimensional Schrodinger equation. They are 
expressed in terms of the matrix coefficients of the separation equations (|l7|)-(|T8l). 

A promising development of the research is classification and study of superintegrable 
(admitting sufficiently many higher symmetries) cases of Pauli equation. Notice that the 
notions of separability and superintegrability are closely related. By now, superintegrable 
physical systems can be regarded as one of the most intensively developed and significant 
fields of mathematical physics. The problem of classifying superintegrable stationary 
Schrodinger equations with scalar potential has been solved by Winternitz with co-workers 
25| and Evans [p6 for space dimensions n = 2 and n = 3 (see also [|J). They have found 



all potentials that allow for separability of the corresponding Schrodinger equation in 
more than one coordinate system. We intend to modify and generalize this approach to 
(1 + 3)-dimensional Pauli equation ([I]). A study of the problem is in progress now and 
will be reported in our future publications. 
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